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Abstract 


Travel-times  of  acoustic  signals  were  measured  between  a  bottom-mounted  source  near 
Oahu  and  four  bottom-mounted  receivers  located  near  Washington,  Oregon,  and  Cali- 
fornia in  1988  and  1989.  This  paper  discusses  the  observed  tidal  signals.  At  three  out 
of  four  receivers,  observed  travel  times  at  Ml  and  52  periods  agree  with  predictions 
from  barotropic  tide  models  to  within  ±30°  in  phase  and  a  factor  of  1.6  in  amplitude. 
The  discrepancy  at  the  fourth  receiver  can  be  removed  by  including  predicted  effects  of 
phase- locked  baroclinic  tides  generated  by  seamounts. 

Our  estimates  of  barotropic  M2  tidal  dissipation  by  seamounts  vary  between  2  x  1016 
and  1  x  1018  erg-s-1.  The  variation  by  two  orders  of  magnitude  is  due  to  uncertainties  in 
the  numbers  and  sizes  of  seamounts.  The  larger  dissipation  (1  x  1018  erg-s"1)  is  the  same 
order  as  previous  estimates  and  amounts  to  4  %  of  the  total  dissipation  at  Ml. 

Thesis  Supervisor:  Dr.  John  L.  Spiesberger 
Woods  Hole  Oceanographic  Institution 


Acknowledgments 

I  owe  many  thanks  to  my  thesis  advisor,  Dr.  John  Spiesberger,  for  his  guidance.  His 
knowledge,  understanding,  and  insight  enabled  the  completion  of  this  thesis. 

I  thank  a  fellow  Naval  Officer,  Paul  Bushong.  He  laid  the  foundations  for  much  of  this 
work  while  studying  at  MIT/WHOI  from  1985  to  1987. 

I  am  indebted  to  the  United  States  Navy  for  giving  me  the  opportunity  to  pursue  my 
graduate  studies  at  MIT  and  WHOI. 

I  thank  Carolyn,  my  wife,  for  her  enduring  love  and  support. 

Finally,  I  thank  the  Lord  Jesus  Christ,  without  whom  I  could  do  nothing. 


Contents 


1  Introduction  8 

2  Data  Acquisition  11 

3  Data  Analysis  13 

3.1  Barotropic  Model  for  Acoustic  travel-time 13 

3.2  Barotropic  Model  Comparisons  with  Observed  Tidal  Signals 14 

3.2.1  Tidal  Signal  Estimation 14 

3.2.2  Comparison  with  Model 15 

3.2.3  Analysis 16 

3.3  Baroclinic  Tide  Prediction 18 

3.3.1  Phase-Stable  Internal  Tides 18 

3.3.2  Influence  on  Travel- Time  at  Receiver  R4 19 

4  Tidal  Dissipation  21 

4.1  Tidal  Dissipation  at  M2  from  Satellite  Observations 21 

4.2  Tidal  Dissipation  by  Seamounts 22 

5  Conclusions  25 
A  Obtaining  Tidal  Currents  from  Elevations  27 
B  Guyot-Generated  Internal  Tides  33 
C    Baroclinic  Tide  Travel-Time  Model  39 


D    Signal  Processing  41 

D.l    The  Transmitted  Signal 41 

D.2    Post-Detection  Processing 42 

D.2.1     Demodulation  and  Averaging 42 

D.2.2    Replica  Correlation 43 

D.2. 3     Phase  Change  Between  Sequence  Records 43 


List  of  Figures 


1-1  Source  and  Receivers 9 

2-1  Available  Tomography  Data 12 

3-1  Phase- Difference  Power-Spectra 15 

3-2  Data  to  Barotropic  Model  Comparisons 17 

3-3  R4  Data-to- Model  Comparisons 20 

A-l  M2  Currents  in  the  NE  Pacific 30 

A-2  JK1  Currents  in  the  NE  Pacific 32 

B-l  Two-Layer  Ocean  Model 34 

B-2  R4  Prediction  Error  as  a  Function  of  Reduced  Gravity 38 


List  of  Tables 


3.1  Barotropic  Model  Signal  Amplitudes 16 

3.2  M2  and  52  Internal  Tide  Corrected  Travel-Times  at  Receiver  R4 19 

3.3  N2  and  K2  Internal  Tide  Corrected  Travel- Times  at  Receiver  R4 20 

4.1  Seamount  Distribution  Models 23 

4.2  Seamount  Dimensions 23 

A.l  Major  Tidal  Frequencies 28 

B.l  Generalized  Two-layer  Internal  Tide  Solution 36 

B.2  Paticular  Two-layer  M2  and  S2  Internal  Tide  Solutions 37 


Chapter  1 
Introduction 


We  measure  travel-times  of  acoustic  signals  between  a  bottom-mounted  source  near  Oahu 
and  four  bottom-mounted  receivers  located  near  Washington,  Oregon,  and  California 
in  1988  and  1989  (Figure  1-1).  The  power-spectra  of  basin-scale  travel-times  exhibit 
prominent  variability  at  tidal  periods  [Spiesberger  et  al.,  1989b].  There  are  two  important 
reasons  for  further  study  of  the  tidal  signals.  First,  numerical  models  for  barotropic 
tides  can  be  validated  by  using  these  new  observations.  Second,  it  is  important  to  detide 
tomography  data  to  better  study  other  ocean  fluctuations  at  similar  periods  [Spiesberger 
et  al.,  1989b]. 

The  tides  can  be  separated  into  barotropic  and  baroclinic  modes,  and  their  associated 
influences  on  travel-time  can  be  separated  into  current  and  displacement  effects.  The  tidal 
currents  affect  travel-time  by  alternating  between  alliance  and  opposition  to  the  direction 
of  sound  propagation.  The  rise  and  fall  of  the  water  column  affects  the  travel-time  by 
causing  periodic  vertical  iso-therm  displacements  and  pressure  fluctuations.  Worcester 
et  al.  [1990]  has  shown  that  barotropic  current  effects  can  be  quantified  by  using  recip- 
rocal transmissions;  one-way  transmissions  are  used  in  this  experiment,  so  we  quantify 
barotropic  effects  by  comparing  measured  data  to  barotropic  tide  model  predictions. 

We  find  that  acoustic  travel-times  computed  from  barotropic  tidal  models  are  similar 
to  observations  at  three  out  of  four  receivers.  The  discrepancy  at  the  fourth  receiver 
can  be  understood  by  including  the  effects  of  phase-locked  baroclinic  tides  generated  by 
seamounts  near  the  transmission  path.    Generation  mechanisms  are  important  both  for 
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Figure  1-1:  The  acoustic  tomography  experiments  of  1988  and  1989  consist  of  a  source 
located  near  Kaneohe  Bay,  Oahu,  and  four  bottom  receivers  (Rl  to  R4)  whose  approximate 
positions  are  shown.  The  distances  between  source  and  receivers  are  all  about  4000  km. 


understanding  the  distribution  of  internal  wave  energy  and  the  energy  budget  of  tidal 
dissipation  [Munk,  1966;  Wunsch,  1975;  Bell,  1975;  Lambeck,  1980;  Hendershott,  1981]. 

This  report  closely  follows  Headrick  et  al.  [1990]  and  is  organized  as  follows.  Chapter 
2  discusses  the  experiment  and  the  estimation  of  acoustic  travel-time.  Comparisons  be- 
tween observations  and  predictions  from  barotropic  and  baroclinic  tide  models  are  made 
in  Chapter  3.  Chapter  4  predicts  the  global  dissipation  of  tidal  energy  at  M2  due  to 
seamounts.  Conclusions  follow  in  Chapter  5. 
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Chapter  2 

Data  Acquisition 


The  travel-time  differences  of  pulse-like  signals  are  measured  at  3  minute  intervals  in 
1988  and  2  minute  intervals  in  1989  between  the  source  at  Oahu  and  the  four  receivers 
(Figure  1-1).  Travel-times  are  intermittently  measured  for  a  total  of  between  four  and 
eight  months  (Figure  2-1).  The  transmission  distances  are  about  4000  km  and  the  transit 
times  are  about  40  minutes. 

Spiesberger  et  al.  [1989a,b]  describes  the  maximum  likelihood  estimator  of  phase 
differences  used  to  estimate  differences  in  travel-time.  The  basic  process  is  outlined  in 
Appendix  D.  With  this  estimator,  the  change  in  acoustic  phase  at  133  Hz,  6<j>,  is  related 
to  the  travel-time  difference,  St,  by, 

>r  =    **  (2.1) 

where  u>c  —  2n  ■  133  (rad-s-1)  is  the  center  radian  frequency  of  the  transmitted  signal,  u  is 
the  radian  frequency  of  the  travel-time  oscillation,  and  At  is  the  geophysical  time  interval 
over  which  travel-time  differences  are  measured  (184  s  in  1988,  122  s  in  1989)  [Spiesberger 
et  al,  1989a]. 
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Figure  2-1:  Shaded  regions  indicate  the  times  when  tomographic  data  were  recorded  at 
each  of  the  receivers. 
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Chapter  3 
Data  Analysis 


3.1      Barotropic  Model  for  Acoustic  travel-time 

The  change  in  acoustic  travel-time  due  to  barotropic  tides  is  given  by, 

Arbt(t,rR)  =  ~[/    Z&t-^-ds  +  h  f    C(s,t-—)ds],  (3.1) 

where  the  geophysical  time  at  the  receiver  is  t,  the  geodesic  between  the  source  and 
receiver  is  Tr,  and  the  mean  sound  speed  is  c0.    The  barotropic  current  and  surface 

displacement  at  location  son  the  geodesic  at  time  /  are  denoted  by  u(s, t ^)  and 

C(s,<  — *p-)  respectively.  The  incremental  on  the  geodesic  is  ds,  and  l{s)  is  the  distance 
from  the  receiver.  The  constant  is  ki  =  3.9  x  10-3  [Spiesberger  et  a/.,  1989b].  The  effects 
of  finite  propagation  time  are  taken  into  account  by  the  correction  term  (  — -^).  The  finite 
propagation  delay  amounts  to  a  phase  correction  of  order  ±10°  at  semidiurnal  periods. 

The  first  term  in  (3.1)  contains  the  effect  of  propagating  sound  with  and  against 
the  tidal  current.  The  second  term,  which  amounts  to  a  1  or  2  ms  correction,  contains 
the  combined  effects  of  vertical  iso-therm  displacement  and  pressure  fluctuations  which 
accompany  the  rise  and  fall  of  the  sea  surface.  For  barotropic  tides  the  travel- time  change 
is  about  the  same  for  different  acoustic  multipaths  [Spiesberger  et  al.,  1989a]. 

The  elevation,  (,  is  obtained  from  the  Naval  Surface  Weapons  Center  (NSWC)  tide 
model  derived  by  Schwiderski  [1978].  Unfortunately,  except  for  M2,  NSWC  models  of  the 
tidal  current  field  are  unavailable.  We  therefore  use  a  linearized  version  of  Schwiderski's 
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[1978]  Laplace  tidal  equations  (LTE)  to  predict  tidal  currents.  The  equations  are, 

%  =  2fi  cos  Bv+j^  [fx(«V  ~  PO]  ~  kdu, 

(3.2) 

|f  =  -20  cos 0u  -  £  [gj(a»7-0O]  "  ^u> 

where  A  and  6  are  longitude  and  colatitude  respectively,  u  and  v  are  east  and  north 

current  speeds  respectively,  77  is  Newton's  equilibrium  tide  (which  is  proportional  to  the 

astronomical  tide  generating  potential),  g  is  the  earth's  gravitational  attraction,  Re  is 

the  radius  of  a  spherical  earth,  0  is  the  Earth's  angular  rate  of  rotation,  and  kj  is  a 

damping  coefficient.  The  constants  a  and  /3  derived  by  Schwiderski  [1978]  correct  for  the 

solid-earth  tide,  ocean  tidal  loading  of  the  solid  earth,  and  the  associated  gravity  potential 

perturbations  of  all  tides  and  loadings.   The  solution  to  (3.2)  for  u  and  v  is  outlined  in 

Appendix  A. 

Equation  (3.1)  is  decomposed  into  contributions  from  each  tidal  frequency,  wjt,  of 

interest,  and  solutions  are  determined  by  collecting  and  integrating  sin  a;*.*  and  cosu;*./ 

terms  separately.  The  resulting  travel-time  change  model  is, 

K 
ATbt(*,rR)  =  £%,t(rR)cosu;fc*  +  ffc(rR)sinu^,  (3.3) 

A:=l 

where  c^Tr)  and  ffc(rR)  are  the  solutions  to  the  separated  integrals,  and  K  is  the  total 
number  of  tidal  frequencies  calculated. 

3.2      Barotropic  Model  Comparisons  with  Observed  Tidal 
Signals 

3.2.1      Tidal  Signal  Estimation 

Power-spectra  of  the  differences  in  acoustic-phase  exhibit  prominent  variability  at  tidal 
periods  [Spiesberger  et  al,  1989a,b]  (Figure  3-1).  For  each  receiver,  tidal  signals  are 
estimated  by  fitting  the  8  most  prominent  tidal  components  to  the  data  (Appendix  A, 
Table  A.l).  Specifically,  we  determine  a^  and  b^;  k  =  1,2, ...  ,8,  in, 

8 

6<tftn,TR)=  ]>>fc(rR)cos^n  +  bfc(rR)sinu^n;  n=  1,2,.. .,N,  (3.4) 
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Figure  3-1:  Power  spectral  density  estimate  for  the  133  Hz  phase-differences  at  receiver 
R2.  This  periodogram  uses  a  Bartlett  data  window  with  42  days  of  data  obtained  in  1988. 
The  tidal  peaks  at  the  A'l,  M2,  and  S2  frequencies  have  rms  values  which  correspond  to 
travel-time  change  amplitudes  of  8.6,  23.5,  and  12.3  milliseconds  respectively. 

by  linear  least  squares  [Lawson  and  Hanson,  1974],  where  the  phase-differences  are  mea- 
sured at  times  tn,  and  the  total  number  of  measurements  equals  N.  Depending  on  the 
receiver,  N  equals  55,000  to  135,000,  so  the  least  squares  problem  is  overdetermined. 
Each  least-squares  fit  to  phase-differences  is  compared  to  barotropic  model  travel-times 
to  determine  how  much  of  the  signal  can  be  attributed  to  barotropic  tidal  effects. 

3.2.2      Comparison  with  Model 

Assuming  all  of  the  (3.4)  signal  is  due  to  barotropic  influences,  each  tidal  component 
is  proportional  to  the  time  derivative  of  the  corresponding  (3.3)  model.  Comparisons 
between  observations  and  the  barotropic  model  are  illustrated  by  the  phasors  of  Figure 
3-2.  Each  phasor  magnitude  corresponds  to  a  data-to-model  signal  amplitude  ratio,  which 
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Rl 

R2 

R3 

R4 

Frequency 

Amp(ms) 

Amp(ms) 

A] 

mp(ms) 

Amp(ms) 

M2 

26 

22 

15 

7 

S2 

11 

7 

4 

3 

K\ 

10 

7 

1 

3 

Table  3.1:  Predicted  acoustic  signal  amplitudes  for  barotropic  tides.  The  model  ampli- 
tudes, (e£  +  fjt)1/2,  are  determined  from  the  constants  in  (3.3)  and  have  rms  errors  of 
about  2  ms. 


from  (2.1)  is, 

2*(4  +  H)1/2 


u;cu,,A*(e2  +  f2)i/2' 
and  each  phasor  angle  corresponds  to  a  difference  between  data  and  model  signal  phases, 

which  is  given  by, 

.      _i  f*k\       .      _i  fek 
tan         : —     —  tan 


3.2.3     Analysis 


b/t  /  V  ffc 


First,  we  discuss  the  comparison  at  two  semidiurnal  periods  (M2  and  S2).  Predicted  tidal 
amplitudes  vary  from  3  to  26  ms  at  our  receivers  (Table  3.1).  At  receivers  Rl,  R2,  and 
R3,  observed  amplitudes  are  0.8  to  1.6  times  predicted  amplitudes  (Figure  3-2).  At  R4 
the  observed  amplitudes  are  2.4  and  3.4  times  predictions.  Predicted  and  observed  phases 
are  within  ±30°  at  Rl,  R2,  and  R3  but  differ  by  up  to  47°  at  R4.  Data  at  Rl,  R2,  and  R3 
agree  with  predictions  from  the  barotropic  model  to  first  order,  but  at  R4,  the  barotropic 
model  is  inadequate.  We  return  to  this  problem  in  the  next  section. 

We  next  discuss  the  comparisons  for  one  semidiurnal  period  {K\).  The  models  for 
barotropic  currents  at  diurnal  periods  are  inaccurate  within  a  ±10°  band  centered  at  30° 
N  (Appendix  A).  The  transmissions  to  all  receivers  intersect  this  band,  so  comparisons 
between  predictions  and  observations  are  problematic.    Nonetheless,  predictions  at  Rl 
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Figure  3-2:  Comparison  of  phase  and  amplitude  of  measured  and  predicted  tidal  signals 
at  four  receivers  (Rl  to  R4)  and  three  frequencies  (M2,  52,  Kl).  Predicted  tidal  signals 
are  obtained  from  the  barotropic  model  given  by  (3.1).  Each  vector  angle  corresponds  to 
a  measurement  phase  that  leads  the  predicted  phase  by  the  indicated  amount.  The  unit 
circle  in  the  upper  left  hand  corner  illustrates  the  vector  angle  convention.  Vector  magni- 
tude corresponds  to  the  ratio  of  measured  amplitude  to  predicted  amplitude.  Predicted 
magnitudes  are  listed  in  Table  3.1.  An  exact  match  between  measurement  and  prediction 
is  represented  by  a  vector  of  angle  zero  and  length  one  (tip  just  touching  the  unit  circle  at 
zero  degrees).  Solid  arrows  correspond  to  measurements  of  the  combined  1988  and  1989 
record,  and  the  dashed  arrows  correspond  to  separate  1988  and  1989  data  measurements. 
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and  R2  resemble  observations  at  tidal  period  Kl  (Figure  3-2).  It  may  be  that  Rl  and 
R2  tidal  signals  are  similar  to  predictions  because  their  transmission  paths  travel  smaller 
distances  within  30°  ±  10°  N  than  paths  to  R3  and  R4  (Figure  1). 

3.3      Baroclinic  Tide  Prediction 

3.3.1      Phase-Stable  Internal  Tides 

Apart  from  generation  regions,  baroclinic  tides  are  not  phase-locked  to  the  ATGF,  and 
their  phases  are  incoherent  at  100  km  distance  (approximately  an  internal  tide  wavelength) 
[Hendershott,1981].  In  generation  regions,  baroclinic  tides  may  be  phase-locked  to  the 
ATGF  [Munk  et  al.,  1981].  Open  ocean  evidence  of  this  phase-lock  is  seen  in  a  phase- 
stable  Ml  baroclinic  tide  over  a  large  guyot  east  of  Hawaii  [Noble  et  al.,  1988].  The 
transmission  path  to  R4  passes  within  a  few  hundred  kilometers  of  several  large  seamounts, 
one  of  which  is  indicated  in  Figure  1-1.  We  hypothesize  that  the  semidiurnal  signals  at 
R4  do  not  agree  with  predictions  from  the  barotropic  model  because  of  the  presence  of 
phase-locked  baroclinic  tides  near  the  seamounts. 

The  transmission  paths  to  Rl,  R2,  and  R3  may  be  near  other  large  seamounts,  but 
because  there  are  no  large  discrepancies  between  the  data  and  barotropic  model  predic- 
tions, it  is  unlikely  that  internal  tides  play  a  zero  order  role  in  determining  tidal  signals 
at  these  receivers. 

Appendix  B  summarizes  a  model  for  the  generation  of  baroclinic  tides  from  the  inter- 
action between  barotropic  currents  and  seamounts.  For  baroclinic  tides,  acoustic  travel- 
times  are  modulated  by  sound-speed  fluctuations  associated  with  vertical  straining  of 
iso-tachs  [Spiesberger  et  al,  1989b].  The  effect  for  the  jth  acoustic  multipath,  rR(j),  is, 

1     /*  1(f)   8r  (f) 

ATbc(/,rR(j))  =  -2   /         Cint(f,<  -  -^)^^rfr,  (3.5) 

where  the  position  along  the  multipath  is  r,  the  vertical  displacement  of  the  internal  tide 
is  Cint?  and  the  vertical  derivative  of  the  potential  sound-speed  is    cpJ ' . 
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M2  S2 

Amplitude(ms)      Phase(°)  Amplitude(ms)      Phase(°) 


Barotropic  Model       7  113  3  324 

Baroclinic  Model      20  117  5  265 


Combined  Model      27  116  7  286 

Measured  Data         24  ±  1 114  ±1  7  ±0.1 279  ±  1 

Table  3.2:  Comparison  of  predicted  and  observed  travel-time  changes  at  M2  and  52  for 
receiver  R4.  The  uncertainties  on  the  measured  data  are  based  on  the  differences  between 
1988  and  1989  least-square  fits.  Phases  are  based  on  a  s\nu>t  time  dependence  with  t  =  0 
at  0000  Greenwich  mean  time  on  1  January  1989. 

3.3.2      Influence  on  Travel-Time  at  Receiver  R4 

Being  the  southern-most  receiver,  R4  has  a  more  east- west  oriented  ray  path  than  the  other 
receivers.  This  orientation  results  in  a  weaker  influence  from  the  generally  north-south 
oriented  M2  and  S2  barotropic  tidal  currents,  so  the  effects  of  phase-locked  baroclinic 
tides  are  more  apparent.  The  predicted  Ml  baroclinic  tide,  generated  from  one  of  several 
seamounts  near  the  transmission  path  to  R4  (Figure  1-1),  has  an  amplitude  and  wave- 
length of  6.2  m  and  98  km  respectively  (Appendix  B)  and  produces  a  21  ms  amplitude 
travel-time  signal  (Appendix  C). 

The  sum  of  the  predicted  M2  baroclinic  tidal  signals,  generated  from  all  prominent 
seamounts  along  the  transmission  path,  is  added  to  the  barotropic  prediction.  The 
barotropic  plus  baroclinic  prediction  for  R4  is  almost  identical  to  observations  at  M2 
and  S2  (Figure  3-3  and  Table  3.2).  The  prediction  is  improved  at  K2  but  not  at  N2 
(Table  3.3).  The  frequency  at  JV2  is  close  to  M2  (Table  A.l),  and  the  observed  amplitude 
at  N2  is  small.  We  suspect  that  the  barotropic  plus  baroclinic  prediction  does  not  fit 
observations  at  AT  2,  because  strong  baroclinic  tides  at  M2  may  be  smeared  into  N2  by 
dynamic  processes  [Hendershott,  1981]. 


If) 


Barotropic 
Model 


M2 


S2 


Combined 
Model 


Figure  3-3:  The  phasors  on  the  left  are  identical  to  the  combined  1988  and  1989  R4  record 
phasors  shown  in  Figure  3-2.  The  phasors  on  the  right  present  the  same  analysis  using 
the  combined  barotropic  and  baroclinic  models  in  Table  3.2. 


N2 

Amplitude(ms)  Phase(°) 

Barotropic  Model      1  37 

Baroclinic  Model      6  315 

Combined  Model    ~6  325 

2  ±  1  331  ±  38 


K2 

Amplitude(ms)  Phase(°) 

~T~  '  153 

3  117 


Measured  Data 


4 

5±2 


126 

153  ±  18 


Table  3.3:  Same  as  Table  3.2,  except  comparisons  are  at  N2  and  A'2  frequencies. 
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Chapter  4 

Tidal  Dissipation 


The  effect  of  the  ocean  tidal  bulge  on  the  Earth-Moon  system  is  a  well  known  phenomena 
that  is  discussed  extensively  by  Lambeck  [1980].  The  friction-induced  phase  lag  between 
£  and  r\  produces  a  torque  in  the  Earth-Moon  system  that  is  responsible  for  the  length- 
ening of  the  day  and  recession  of  the  Moon  at  rates  of  about  1.8  ms-century-1  and  3.7 
m-century-1  respectively  [Christodoulidis,  et  al.  1988].  The  applied  torque  transfers  an 
amount  of  orbital  momentum  to  the  Moon  equal  to  the  rotational  momentum  lost  by  the 
Earth  [Lambeck,  1977].  The  net  result  of  this  momentum  transfer  is  an  overall  reduction 
in  system  kinetic  and  potential  energy.  System  energy  is  dissipated  for  the  most  part 
in  oceanic  tides,  with  some  fraction  (which  has  yet  to  be  quantified)  occurring  through 
conversion  to  baroclinic  tides  [Lambeck,  1980].  In  this  chapter  we  estimate  the  barotropic- 
to-baroclinic  conversion  rate  of  the  M2  tide  at  seamounts. 

4.1      Tidal  Dissipation  at  M2  from  Satellite  Observations 

The  rate  of  energy  dissipation  is, 

where  —^-  is  the  total  tidal  acceleration  of  the  Earth  and  ^  is  the  acceleration  of  the 
Moon's  orbit  in  rad  •  s~     [Lambeck,  1977].    Recent  estimates  for  Ml  accelerations  have 
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been  obtained  from  the  orbital  perturbations  of  17  satellites,  and  have  values  given  by, 

(4.2) 


,ni2 


[Christodoulidis,  et  al.,  1988].  Substituting  (4.2)  into  (4.1)  yields  2.5  x  1019  erg-s-1  at  M2 
(about  75  %  of  all  tidal  dissipation). 

4.2      Tidal  Dissipation  by  Seamounts 

We  consider  two  models  for  the  distributions  of  seamounts.  In  the  first  model  the  number 
of  seamounts  with  a  basal  radius  r  >  R  (km)  is  given  by, 

N(R)  =  1.44  x  106exp[-0.63£].  (4.3) 

Equation  (4.3)  is  an  extrapolation  of  Jordan  et  al.  's  [1983]  results,  for  the  eastern  Pacific, 
to  the  world  ocean.  In  the  second  model,  the  number  of  seamounts  is  estimated  from 
altimeter  measurements  by  Seasat  [Craig  and  Sandwell,  1988].  Seasat  detects  many  large 
seamounts  which  are  absent  in  the  model  by  Jordan  et  al.  [1983],  but  it  does  not  detect 
small  seamounts,  as  the  altimeter  profile  sensitivity  to  seamounts  falls  off  rapidly  below  a 
basal  radius  of  20  km  [Craig  et  al.,  1988].  These  two  seamount  distributions  are  outlined 
in  Table  4.1. 

The  rate  of  energy  conversion  from  barotropic  to  baroclinic  tides  at  each  seamount  is 
approximately, 

^  =  2(Ax)cgiYl,  (4.4) 

where  Ax  is  the  horizontal  dimension  of  each  seamount  perpendicular  to  the  direction  of 
barotropic  current,  cg{  is  the  group  speed  of  the  internal  tide  (Equation  (B.22)),  and  E{  is 
the  average  energy  density  of  the  internal  tide  (Equation  (B.21)).  The  dimension,  Ax,  is 
obtained  from  the  seamount 's  basal  radius  according  to  Table  4.2.  If  we  take  the  barotropic 
current  speed  to  equal  2  cm-s-1  and  the  inertial  frequency  to  equal  0.73  X  10~4  s_1  (30° 
N),  then  equations  (4.4)  and  (B.23)  yield  a  total  energy  conversion  rate  of  2  x  1016  and 
1  X  1018  erg-s-1  for  the  exponential  and  Seasat  distributions  respectively.  These  estimates 
are  proportional  to  Ax  and  the  squares  of  both  seamount  height  and  the  barotropic 
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Number  of  Seamounts 

Basal  Radius 

Exponential 

Seasat 

r(km) 

Model 

Model 

1-2 

360,000 

0 

2-4 

290,000 

0 

4-6 

83, 000 

0 

6-10 

30,000 

200 

10-  15 

2500 

1200 

15-20 

110 

6000 

20-26 

5 

8800 

26-32 

0 

6800 

32-38 

0 

4800 

38-44 

0 

2800 

44-50 

0 

1600 

Table  4.1:  Seamount  distribution  for  the  world's  oceans  as  a  function  of  basal  radius  by 
Jordan  et  al.  [1983]  (exponential  model)  and  Craig  and  Sandwell  [1988]  (Seasat  model). 
The  Seasat  model  is  obtained  by  converting  Craig  and  SandwelVs  [1988]  histogram  of 
peak-to-trough  separation,  p,  to  basal  radius  using  r=(0.83)p  in  accordance  with  mea- 
sured seamount  slopes.  The  numbers  so  obtained  are  multiplied  by  four  to  account  for 
undetected  seamounts  [Craig  et  al,  1988].  Some  of  the  discrepancy  between  models  may 
be  due  to  large  seamounts  in  the  western  Pacific  that  were  not  represented  in  the  Jordan 
et  al.  [1983]  survey. 


Rectangular 

Seamount  Dimensions 

Dimension 

Exponential 
Model 

Seasat 
Model 

Ax 

Az 

(0.62)r 
(0.21)r 

(0.62)r 
2500  m 

Table  4.2:  Seamount  dimensions  for  the  exponential  and  Seasat  models  (Table  4.1).  The 
length  and  width  of  each  seamount  is  Arc,  and  the  height  is  Az.  Basal  radius  is  denoted 
by  r.  The  scale  factor  for  Ax,  (0.62),  is  derived  from  a  mean  flatness  of  0.31  reported  for 
85  seamounts  in  the  Pacific  [Smith,  1988].  The  Az  scale  factor  for  the  exponential  model 
is  derived  from  the  mean  height-to-basal-radius  ratio  of  0.21  [Smith,  1988].  The  height  is 
set  to  2500  m  for  the  Seasat  model. 
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current  speed  (Equation  (B.23)).  The  exponential  model  yields  insignificant  dissipation 
because  it  contains  relatively  few  wide  seamounts  (Table  4.1).  The  Seasat  model  yields  a 
dissipation  that  is  about  4  %  of  the  total  at  M2.  Munk  [1966]  and  Bell  [1975]  obtained 
rates  of  5.0  X  1018  and  2.5  X  1018  erg-s-1  respectively  for  the  conversion  to  baroclinic  tides 
by  bottom  roughness.  These  estimates  are  the  same  order  as  we  obtain  using  the  Seasat 
distribution. 
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Chapter  5 
Conclusions 

We  reach  the  following  conclusions: 

1.  Because  barotropic  signal  predictions  at  Rl,  R2,  and  R3  compare  favorably  with 
observations  and  because  barotropic  wave-lengths  have  basin-scales,  the  principle 
discrepancy  between  barotropic  signal  predictions  and  observations  at  R4  is  probably 
not  due  to  errors  in  Schwiderski's  [1978]  model.  Furthermore,  because  barotropic 
plus  baroclinic  signal  predictions  compare  favorably  with  observations  at  R4,  we 
believe  the  Schwiderski  [1978]  model  is  adequate  for  barotropic  predictions. 

2.  Seamounts  may  be  responsible  for  a  small  (~  5  %),  but  non-negligible,  fraction  of 
the  dissipation  of  barotropic  tides.  The  dissipation  is  sensitive  to  the  heights,  radii, 
and  numbers  of  seamounts  (all  poorly  known). 

3.  The  generation  of  baroclinic  tides  at  seamounts  can  be  accounted  for  with  a  two- 
layer  model  of  the  ocean  (Appendix  B). 

4.  Barotropic  tidal  currents  can  be  predicted  from  Schwiderski's  [1978]  tidal  elevation 
predictions  with  our  model  (Appendix  A),  except  within  ±10°  of  30°  North  and 
South  (diurnal  frequencies)  and  near  90°  N  (semidiurnal  frequencies). 

5.  Comparison  of  measured  tidal  signals  (from  basin-scale  acoustic  transmissions)  with 
predicted  tidal  signals  (from  barotropic  models)  can  be  used  to  detect  baroclinic 
tides  which  are  phase-locked  to  the  ATGF. 
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6.  It  is  difficult  to  make  a-priori  predictions  of  tidal  signals  in  basin-scale  transmissions, 
because  generation  regions  of  baroclinic  tides  are  poorly  mapped. 
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Appendix  A 

Obtaining  Tidal  Currents  from 
Elevations 


We  solve  (3.2)  for  u  and  v  by  defining, 


u  —  b\  cosujt  +  62  sin u>t, 
v  —  a\  cosu;£  +  a2sinu;£, 


and. 


(A.l) 


(A.2) 


jdne  [w\(a7l  ~  PO]  =  yi  cosujt  +  2/2  smut, 

-  ji  [§§(QV  -  P0\  =  *\  cos ut  +  z2  sin  wt, 
where  u>  is  the  tidal  frequency  of  interest  and  ai,  a2,  61,  62,  t/i,  t/21  z\,  and  z2  are  constants. 
The  elevation,  £,  and  Newton's  equilibrium  tide,  77,  are  periodic  with  frequencies  shown 
in  Table  A.l.  The  j/'s  and  z's  are  computed  from  (A.2)  by  expanding  known  functions,  £ 
and  77,  in  terms  of  sines  and  cosines.  Substituting  (A.l)  and  (A.2)  in  (3.2)  leads  to  four 
linear  equations,  in  the  unknowns  ai,  <*2,  61,  and  62,  which  can  be  written  in  matrix  form 


0 

-2ft  cos 

0 

—u> 

kd 

—  U) 

kd 

0 

2Qcos 

2Q  cos  0 

0 

kd 

uj 

kd 

u 

2ft  cos 

6 

0 

wa,  \ 


('2 
V    fc2    / 


'^ 


*2 


(A.3) 
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Semidiurnal  Modes 

Diurnal  Modes 

Frequency 

Frequency 

Component 

(cpd) 

Component 

(cpd) 

M2 

1.93227 

Kl 

1.00273 

S2 

2.00000 

01 

0.92954 

N2 

1.89600 

PI 

0.99726 

K2 

2.00547 

Ql 

0.89325 

Table  A.l:  Major  Tidal  Frequencies 


The  solution  to  this  matrix  equation  determines  u  and  v. 
We  use  the  values, 

a  =  .69,    (3  =  .90,  (A.4) 

derived  by  Schwiderski  [1978],  and  the  value  of  the  linear  damping  coefficient,  k^,  is 
selected  by  minimizing  the  phase  difference  between  predictions  of  AT(^,M2,r^3)  using 
Schwiderski  [1983]  and  damped  LTE  currents.  The  resulting  best  value  is, 

kd  =  9.0x  10-V1.  (A.5) 

The  value  of  kj  also  controls  the  Q  of  the  LTE  model.  We  estimate  Q  by  comparing 
the  equations  of  motion  to  a  forced  linear-damped  oscillator, 

ft     T  ftT 

-^■  + k—  +  u>£x  =  Acosut,  (A.6) 


Q  =  £,  (A.7) 


whose  Q  is  defined  as 

k 
where  the  resonant  frequency  is, 

"R  =  K2-y)1/2,  (A.8) 

[Marion,  1970].  By  letting  east  be  the  x-direction,  and  substituting  (A. 2)  into  (3.2),  we 

obtain  the  following  x-direction  momentum  equation, 

d  x  dx 

-Qp  +kd-7?-- fv  =  yi  cos ut  +  y2 smut,  (A.9) 
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where  /  =  20.  cos  6  is  the  inertia!  frequency.  Since, 

Iff  |  «1.4xl0-6, 

\fv\  «  1.4  x  10~6, 

(»?  +  Vlf'2  ~4.4xl0-7, 

\kdu\  w  9.0  x  10~8, 

the  dominant  terms  in  (A.9)  can  be  written 

du 


(A.10) 


gt      f»  =  0.  (A.11) 


Substituting  the  solution  to  (A. 11), 

i 
and  the  relation, 

into  (A9),  we  obtain  the  approximation, 
82x 


ILOU 

T' 


(A.12) 


(A.13) 


,  kd—  +  u  x  ~  t/i  cosojt  +  y2  smut.  (A. 14) 


^-^^_U1X10-V' 
^~  t.  -I.."-    onvm-6c-i    _  i0'  ^A'15^ 


Comparing  (A.6,7,8)  with  (A.  14)  we  see, 

2 

fcd  ~  Ard        9.0  X  10"6f 

where  1.41  x  10_4s_1  is  the  frequency  of  the  M2  tide.  This  agrees  with  estimates  ranging 
from  5  to  25  provided  by  Hendershott  [1981]. 

A  comparison  of  the  damped  LTE  model  with  the  Schwiderski  [1983]  model  for  Ml 
in  the  North  Pacific  ocean  demonstrates  good  agreement  (Figure  A-l).  Both  models  also 
generally  agree  with  the  data  presented  by  Luyten  and  Stommel  [1990]  extracted  from 
deep  moored  current  meter  records.  By  extension,  we  believe  that  the  damped  LTE 
model  should  be  an  accurate  first  order  predictor  of  the  52,  A^2,  and  A'2  semidiurnal 
currents. 

The  damped  LTE  model  does  not  perform  as  well  as  the  Schwiderski  [1983]  Ml  current 
model  in  the  vicinity  of  amphidromic  points.  The  finite  difference  approximations  we  use 
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Figure  A-l:  The  M2  barotropic  tidal  current  in  the  northeast  Pacific  Ocean  at  0600 
Greenwich  mean  time  on  1  January  1989  as  determined  from  (a)  the  Schwiderski  [1983] 
model  (solid  arrows)  (b)  linear-damped  Laplace  tidal  equations  (Equation  (3.2);  dashed 
arrows)  (c)  data  extracted  from  deep  moored  current  meter  records  (two  arrows  desig- 
nated by  circles  around  their  tails).  The  current  meter  analysis  is  taken  from  Luyten  and 
Stommel  [1990] 
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for  the  spatial  derivatives  in  (3.2)  are  not  accurate  near  these  areas  of  rapidly  changing 
phase.  These  inaccuracies  are  limited  to  a  100  km  (the  grid  size  of  Schwiderski's  model) 
radius,  and  the  net  effect  on  predicted  travel-time  change  is  less  than  5  %  in  amplitude 
and  3  degrees  in  phase  for  a  4000  km  geodesic  passing  through  an  amphidrome. 

The  damped  LTE  model  has  an  additional  difficulty  for  diurnal  currents  within  a  ±10° 
band  centered  at  the  inertial  latitude, 

0}  =  sin-1(^/2ft).  (A. 16) 

For  diurnal  tides,  u;^  ~  1  cpd  and  Of  =  30°  N,  and  for  semidiurnal  tides,  u>k  ~  2  cpd  and 
Of  =  90°  N.  The  difficulty  can  be  seen  by  setting  kj  =  0  and  solving  (3.2)  for  u  and  v. 
The  resulting  equations, 

g  [-2fi  cost  jL(ar,-K)+^  £(ai>-/3<)] 
U  ~  i?(4Q2cos20-u,2)  ' 

(A.17) 


R(4U2  cos2  6-w2)  ' 

have  denominators  equal  to  zero  at  the  inertial  latitude.  These  zeros  are  canceled  by 
numerator  zeros  when  u,  v,  and  (  are  determined  exactly  from  the  same  equations  of  mo- 
tion, but  a  singularity  occurs  when  the  numerators  are  approximated  by  finite  differences. 
The  finite  kd  we  use  in  (3.2)  removes  the  singularity  in  (A.17),  but  the  solution  is  still 
inaccurate  near  30°  N  (Figure  A-2). 
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Figure  A-2:  The  K\  barotropic  tidal  current  in  the  northeast  Pacific  Ocean  at  0600 
Greenwich  mean  time  on  1  January  1989  as  determined  from  the  linear-damped  Laplace 
tidal  equations  (Equation  (A.2)).  Note  the  irregular  performance  of  the  model  in  the 
vicinity  of  30.1°  N. 
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Appendix  B 

Guyot- Generated  Internal  Tides 


In  order  to  determine  the  coupling  between  barotropic  tidal  currents  and  internal  tides 
generated  by  guyots,  we  use  the  two-dimensional  two-layer  model  illustrated  by  Figure  B- 
1.  This  model  is  applicable  to  sea  plateaus  whose  major  horizontal  extent  is  perpendicular 
to  the  direction  of  the  barotropic  current.  We  follow  the  procedure  used  in  deriving 
the  coupling  between  surface  tides  and  coastal  generated  internal  tides  [Rattray,  I960]. 
Variables  are  a  function  only  of  the  distance  from  the  center  of  the  plateau  (y-direction) 
in  this  two-dimensional  model. 

A  rigid-lid  approximation  is  used.  The  momentum  equations  are, 

dtv'  +  fu'  =  0,  (B.l) 

8tv"  +  fu"  =  -grdy(",  (B.2) 

and, 

8tu'  -  fv'  =  0, 

„  „  <B'3) 

dtu    -  fv"  =  0, 

where  u  is  the  current  parallel  to  the  plateau,  (')  terms  correspond  to  upper  layer  and  (") 
to  lower  layer  terms,  /  is  the  inertial  frequency,  and, 

gr  =  gAp/p,  (B.4) 

is  reduced  gravity,  where, 

Ap  =  p"-P'.  (B.5) 
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Figure  B-l:  Definition  Figure  for  a  two-layer  ocean  model  used  to  predict  the  coupling 
between  a  barotropic  tidal  current  and  the  internal  tide  generated  by  a  seamount  whose 
extent  perpendicular  to  the  direction  of  the  current  is  considered  infinite. 


Subtracting  (B.2)  from  (B.l),  we  find, 

dt(v'  -  v")  +  f{u'  -  u")  =  grdy(", 
and  the  continuity  equations  for  each  layer  are, 

dyQl'v')   =   dtC, 

dy(h"v")  =  -dt(". 


(B.6) 


(B.7) 


Let  vs  =  (a/d)expi(u>t  +  </>),  which  is  essentially  constant  for  any  distance  y  small  com- 
pared to  the  wavelength  of  the  surface  tide.  We  have  contributions  to  the  currents  in  the 
two  layers  from  both  the  barotropic  current  and  the  currents  associated  with  the  internal 
tide.  Thus, 

v'  —  vs  +  v't       u'  =  us  -f  u\ 

v"  =  vs  +  v"     u"  =  us  +  u", 
and  from  the  rigid-lid  approximation, 


(B.8) 


C'  =  o, 


(B.9) 
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(B.10) 


We  choose  the  following  solutions  which  satisfy  (B.3)  and  (B.7), 

y  < -l  -I  <y  <  I  y>£ 

Ct     -Betk™  Asmk-iy  Be~lk™ 

vs       l/d2  1/d!  l/d2 

v'i       -(u/k2h')Belk?y  -{iw  J kxh') A  cos  k^y  -(u/k2h')Be-tk™ 

v'/      (u>/k2h%)Beik*»  (iu/ kxh'O A  cos  kiy  (cj/k2h'i)Be-lk™ 

us       —ijjud2  -if/ojdi  -if/uid2 

u'i      {iflk2h')Beik™  -Ulkxh')Acoskxy  (if/k2h')Be-tk™ 

<      -(if/k2h%)Beik™  {f/k2h'{) A  cos  kxy  -(if/k2h%)Be-lk™, 
with  the  common  multiplier, 

aexp,'(wH-*).  (B.ll) 

To  satisfy  (B.6), 

*?  =  $?#*>!    J  =  l>2'  (B-12) 

In  order  to  patch  the  solutions  at  the  edges  of  the  plateau,  (,"  must  be  continuous  at  y  =  I 
and  y  =  -£.  Thus, 

Be~ik2e  =  AsmkxL  (B.13) 

The  solutions  (B.10)  already  satisfy  the  conservation  of  total  mass  flux  h'v'  -f  h"v"  = 
dvs  +  h'v'i  +  h"v"  across  y  =  I  and  y  =  —£,  but  conservation  of  p'  mass  flux  requires  that 
h'v'  =  h'(vs  +  v'j)  be  continuous  also.  To  satisfy  this  requirement  we  let, 

A  =  Aoe-^~^\ 
B  =  B0e-i(a-*f2-k*e). 
Substituting  (B.14)  into  (B.13)  leads  to, 

A0smk1i=  B0.  (B.15) 

Using  these  relations  and  conservation  of  p'  mass  flux  we  find, 

h'       u>  h'       iuj 

— Ao(sina  +  i  cos  a)  sin  k\t  —  - —  Ao(sina  +  icoso)cos/:i£.  (B.16) 

d2       k2  d\       k\ 

By  equating  the  imaginary  part  of  (B.16)  we  obtain, 

tano^-p-  tan  k\l,  (B.17) 

k2 
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(B.14) 


Sint 

aAosm  k\y 

<f>  -  a  +  90 

vs 

a/d\ 

4> 

< 

WjjAq  cos  kiy/kih' 

4>  —  a 

< 

oujAq  cos  k\  y/ki  h'{ 

<f>  -  a  +  180 

Us 

af/codi 

0  +  270 

< 

af  Aq  cos  kiy/k^h' 

0  -  a  +  270 

< 

afAo  cos  kiy/kih" 

0  -  a  +  90 

|y|<* M>* 

Amplitude  Phase(°)  Amplitude Phase(°) 

aB0(\y\/y)  <f>-a  +  90-k2(\y\-e) 

a/d2  <t> 

awB0/k2h'  <j>  -  a  +  270  -  k2(\y\  -  I) 

au>B0/k2h'2'  4>  -  a  +  90  -  k2(\y\  -  I) 

af/u>d2  <f>  +  270 

afB0/k2ti  <j>  -  a  +  180  -  k2(\y\  -  i) 

afB0/k2h2'  4>-  a-  k2(\y\  -  I) 

Table  B.l:  Generalized  two-layer  internal  wave  solution  for  a  barotropic  north-south  tidal 
current  propagating  over  an  undersea  plateau  of  infinite  east-west  extent. 


and  by  equating  the  real  part,  we  can  write, 

.  1         1  ,  k\h'  cosa  .„    ^N 

A0  =  (-, r)~ r^r,  B.18 

yd2      d/wcosM  v         ; 

which  leads  to, 

Bo  =  (^--h—  sma.  (B.19) 

d2       d\      u> 

The  result  of  this  solution  is  a  standing  internal  wave  over  the  extent  of  the  plateau  which 
patches  to  traveling  internal  waves  at  either  end.  The  generalized  solution  is  collected  in 
Table  B.l,  and  the  particular  solution  for  aguyot  located  at  30.3°N,  222.5°E  is  summarized 
in  Table  B.2. 

The  cycle-average  baroclinic  energy  density  for  this  two-layer  ocean  is, 


t70    2pl{Ch 


Ygr  +  h'(utf  +  h'(vtf  +  ti\u''Y  +  h"(v'/Y]dt,  (B.20) 


where  T  is  the  period  of  the  internal  tide.    Substituting  the  Table  B.l  amplitudes  and 
phases  in  (B.20)  results  in, 

El  =  \p9r(vsd2B0)2  [jjrrp) >  (B-21) 

where  the  resulting  energy  density  is  in  joules  per  square  meter.  This  baroclinic  energy 
propagates  at  a  group  speed  that  from  (B.12)  is  given  by, 

u>2  -f2 

(B.22) 


k2u 
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M2 

52 

Greenwich 

Greenwich 

Constituient 

Amplitude     Phase(°) 

Amplitud< 

>     Phase(°) 

6.21  m 

166  +  k2y 

2.08  m 

202  +  k2y 

Vs 

2.3  cm-s- 

68 

0.7  cm-s" 

104 

< 

1.9  cm-s- 

346  +  k2y 

0.6  cm-s-: 

22  +  k2y 

V» 

0.3  cm-s- 

1      166  +  k2y 

0.1  cm-s-] 

202  +  k2y 

us 

1.2  cm-s- 

1      158 

0.4  cm-s" 

194 

< 

1.0  cm-s- 

1     76  +  k2y 

0.3  cm-s- 

112  +  k2y 

u'! 

0.2  cm-s- 

1      256  +  k2y 

0.1  cni'S" 

292  +  k2y 

Table  B.2:  Two-layer  Ml  and  52  internal  tide  solutions  for  regions  north  of  a  22  by  111 
km  undersea  plateau.  Seamount  is  located  at  30.3°N,  222. 5°E  (Figure  1-1),  I  =  11  km, 
dx  =  2200  m,  d2  =  4700  m,  h'  =  700  m,  gr  =  .0059  m-s"2,  k2(M2)  =  6.41xl0~5  m"1 
(wavelength  is  98  km),  and  k2(S2)  =  6.71xl0-5  m_1  (wavelength  is  94  km). 


The  group  speed,  energy  density,  and  horizontal  extent  of  each  seamount  are  combined 
in  (4.4)  to  obtain  the  barotropic-to-baroclinic  tidal  energy  conversion  rate.  Expansion  of 
(4.4)  in  terms  of  fundamental  parameters  results  in, 

d£  A.  /  ^-N2/,,2..f2\1/2 

(It 


>o^-4-i)2r- 


p 


(B.23) 


where. 


tana 


dxh'^11' 
d2h'[ 


tan 


f'2 


grh'h\ 


1/2 


(B.24) 


The  internal  tide  in  Table  B.2  has  an  energy  density  and  group  speed  of  about  1.5  X  105 
erg  •  cm-2  and  160  cm-s-1  respectively.  The  energy  density  near  this  seamount  is  about 
17  %  of  that  contained  in  the  parent  barotropic  tide,  and  the  baroclinic  conversion  rate 
is  about  5.3  X  1014  erg-s-1  (53  Mw). 

The  reduced  gravity  (B.4)  used  in  the  dispersion  relation  (B.12)  is  subject  to  the 
choices  for  the  representative  densities  p'  and  p"  in  the  upper  and  lower  layers.  Values  of 
gr  of  order  10-2  m-s-2  are  typical  for  two-layer  ocean  models,  but  using  gr  =  10-2  m-s-2 
is  not  optimal  for  minimizing  the  mean  square  error  between  the  combined  model  and  the 
measured  travel  time  signals  at  R4.  The  optimal  value  of  reduced  gravity  chosen  for  this 
model,  as  illustrated  by  Figure  B2,  is  about  5.9  X  10-3  m-s-2. 
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0.004  0.006         ,       0.008  0.010 

Reduced    Gravity   (m/s2) 

Figure  B-2:   Magnitude  of  Receiver  4  barotropic  plus  baroclinic  signal  prediction  errors. 
The  errors  are  plotted  as  a  function  of  reduced  gravity. 
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Appendix  C 

Baroclinic  Tide  Travel-Time 
Model 


Equation  (3.5)  predicts  the  travel-time  modulation  of  a  single  multipath  caused  by  one 
seamount.  The  multipaths  oscillate  about  the  sound  axis,  and  a  typical  wavelength  (upper 
plus  lower  loop)  is  50  km.  The  potential  sound  speed  gradient  in  the  region  of  the  Table 
B.2  plateau  is  approximately, 

dcp        I    0.13  s-1     above  axis 
0z        I    0.0  s-1       below  axis, 

where  positive  z  is  upward.  We  assume  the  average  ray  path  is  above  axis  for  50  %  of  its 
travel  adjacent  to  the  111  km  length  of  the  plateau,  so  the  average  change  in  travel- time 
at  the  Ml  period  is, 

AibcW-  (2)(1479ms-i)2  ~lle  mS'         (C'2) 

where  0^int  is  the  phase  of  the  M2  internal  tide,  and  u>  is  the  M2  frequency. 

An  evaluation  of  Ns  seamount-generated  tides  along  the  geodesic  between  source  and 
receiver  will  produce  A^  baroclinic  signals  at  each  tidal  frequency.  The  vector  sum  of 
baroclinic  models  plus  the  barotropic  model  results  in  the  combined  model, 

Ns 

ATtot(t,ojk)  =  ATht(t,cjk)  +  £  ATbc(n,*,u;fc),  (C.3) 

n-l 
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where  ATbc(n,2,u;fc)  is  the  baroclinic  signal  prediction  for  the  nth  seamount  at  frequency 
Wfc. 
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Appendix  D 
Signal  Processing 

D.l      The  Transmitted  Signal 

Pulse  compression  techniques  described  by  Spiesberger  et  al.  [1980]  and  Metzger  [1983]  are 
utilized  to  increase  signal-to-noise  ratio  without  degrading  the  phase  difference  resolution. 
The  transmitted  signal  has  a  center  frequency  of  133  Hz  and  a  bandwidth  of  about  17  Hz. 
The  center  frequency  is  phase  modulated  every  8  cycles  using  a  511  digit  maximal  shift 
register  sequence  having  a  period  of  511(8/133)  «  30.7  s  [Spiesberger  et  al,  1989a].  The 
ideal  transmitted  waveform  is, 

s(t)  =  s0cos[2irfct  -  |(rA(t)  -  -)]  (D.l) 

[Spiesberger  et  al,  1980],  where  fc  is  the  center  frequency,  r*.  represents  the  shift  register 
sequence  of  ones  and  zeros,  and, 

OO  £, 

fc(t)  =  ^u(t-—),  (D.2) 

/=o  U6 

where  u(t)  is  the  unit  step  function. 
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D.2      Post-Detection  Processing 

D.2.1      Demodulation  and  Averaging 

The  received  signal,  z(t ),  can  be  approximated  by, 

Z(t)=    £&„*(<  -Tn),  (D.3) 

n=l 

which  is  the  weighted  sum  of  Nmp  multipaths  each  having  gain  bn  and  delay  rn  [Spiesberger 
et  al.,  1980].  Demodulation  of  the  received  signal  is  achieved  by  sampling  at  4  •  133  Hz  = 
532  Hz.  Samples  of  an  ideal  single-path  arrival  are, 

*[n]  =  cos(^)cos[^(r,[n]  -  1/2)]  +  sin(^)sin[|(r%]  -  1/2)],  (D.4) 

where, 

oo 

fc[n]  =  ^u[n-32/],  (D.5) 

/=o 
and  u[n]  is  the  unit  step  sequence.    Because  the  sampling  frequency  is  four  times  the 

carrier  frequency,  the  even  numbered  samples  are  the  cosine  of  the  demodulated  phase 

angle  (weighted  by  ±1)  and  the  odd  numbered  samples  are  the  sine  of  the  demodulated 

phase  angle  (weighted  by  ±1).    The  demodulated  phase  angles  are  low-pass  filtered  by 

block  averaging  to  achieve  a  sequence  of  complex  demodulates  defined  as, 

c[j]  =  x[j]  +  iy[j};  j  =  0, 1, . . . ,  J  -  1,  (D.6) 

where  the  square  root  of  minus  one  equals  i, 

x[j}=  £(-ir+M2(m  +  4j)],  (D.7) 

m=0 

and, 

3 

y[j}=  £(-ir+i*[2(m  +  4j)  +  l],  (D.8) 

m=0 

[Spiesberger  et  al,  1980].  The  resulting  sample  rate  is  four  complex  demodulates  per 
digit,  so  there  are  J  =  4  •  511  =  2044  complex  demodulates  per  sequence.  The  complex 
demodulates  are  averaged  over  6  consecutive  sequences  in  1988  resulting  in  one  averaged 
sequence  record  every  184  seconds.  The  demodulates  are  averaged  over  4  sequences  in 
1989.  The  signal-to-noise  ratio  increases  by  101og106  =  8  dB  following  the  six-period 
average  in  1988  and  by  101og104  =  6  dB  following  the  four-period  average  in  1989. 
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D.2.2      Replica  Correlation 

One  period  of  a  66.5  Hz  sampling  of  (D.l)  phase  angles  (expressed  as  complex  demodu- 
lates) is  contained  in  the  following  sequence  of  length  J, 

*M  =  cos[|(rfcli]  -  \))  +  » "m[|(r*M  -  5)];  J  =  0, 1, ...,/-  1,  (D.9) 

where, 

00 

H?]  =  £u[i-4/],  (D.10) 

and  u[j]  is  the  unit  step  sequence.  The  received  complex  demodulates  (D.6)  are  circularly 
convolved  with  a  sequence- reversed  (matched  filter)  replica  of  the  transmitted  demodu- 
lates (D.9).  The  pulse-compressed  output  is  given  by, 

j-\ 
d[m)  =  £  c[j]s;[j  +  m];  m  =  0, 1, . . .  ,J  -  1,  (D.ll) 

i=o 

where  the  periodic  sequence,  s^j],  is  given  by, 

00 
h[J)=    £    «*Ii  + rJ],  (D.12) 

r=— 00 

[Oppenheim  and  Schafer,  1989],  and  s£  is  the  complex  conjugate  of  s^.  The  replica 
correlation  (D.ll)  compresses  30.7  s  of  energy  along  each  multipath  into  a  single  pulse 
of  one  digit  length  (~  0.06  s)  [Spiesberger  et  al.,  1989a],  whose  amplitude  and  location 
depend  on  the  gain,  6n,  and  delay,  rn,  in  (D.3).  The  signal-to-noise  ratio  increases  by 
101og10511  =  27  dB  following  pulse  compression,  and  a  sliding  average  of  4  complex 
samples  (separated  by  0.015  s)  increases  the  signal-to-noise  ratio  by  another  6  dB  at  the 
cost  of  some  resolution  [Spiesberger  et  al.,  1989a]. 

D.2.3      Phase  Change  Between  Sequence  Records 

The  change  in  phase  of  the  mth  pulse-compressed  demodulate  from  record  r  —  1  to  r  is 
given  by  the  phase  of, 

d[m,r  —  l]d*[m,  r]. 

The  average  change  in  phase  from  record  r  -  1  to  r  is  given  by  the  phase  of 

M 

Aei64>[r]  =  J2  d[m,r-  l]d*[m,r],  (D.13) 
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where  there  are  M  demodulates  containing  a  signal  [Spiesberger  et  ah,  1989a].  The  phase 
in  (D.13)  is  referred  to  as  the  cross-correlated  phase,  because  it  is  the  phase  of  the  complex 
cross  correlation  of  adjacent  pulse-compressed  records  at  zero-lag.  The  cross-correlated 
phase  can  be  converted  to  the  average  travel-time  difference  between  records  using  (2.1). 
The  precision  of  the  resulting  estimation  of  travel-time  difference  is  about  135  /is  [Spies- 
berger et  al,  1989a]. 
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